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Abstract 



We study the Einstein-scalar field system with positive cosmological constant and spherically 
symmetric characteristic initial data given on a truncated null cone. We prove well-posedness, global 
existence and exponential decay in (Bondi) time, for small data. From this, it follows that initial data 
' close enough to de Sitter data evolves to a causally geodesically complete spacetime (with boundary), 

which approaches a region of de Sitter asymptotically at an exponential rate; this is a non-linear 
stability result for de Sitter within the class under consideration, as well as a realization of the cosmic 
no-hair conjecture. 

cr 

H ■ 1 Introduction 

OJj' 

The introduction of a positive cosmological constant A into the Einstein field equations allows one to 
, model inflation periods (large A) as well as the "recent" period of accelerated expansion (small A), and 

consequently plays a central role in modern cosmology. This adds to the relevance of studying initial value 
problems for the Einstein-matter field equations with positive cosmological constant. For such problems 
a general framework is provided by the cosmic no-hair conjecture, which states that generic expanding 
solutions of Einstein's field equations with a positive cosmological constant approach the de Sitter so- 
i ■ lution asymptotically. This conjecture as been proved for a variety of matter models and/or symmetry 

conditions [FrTMl IWal83l IRen04l ITNR03I ITNN051 IRinOSl iRSfM IBey09c[ |Spell| , but the complexity of 
£S| ■ the issue makes a general result unattainable in the near futurej^j 

Here we will consider the spherically symmetric Einstein-scalar field system with positive cosmological 
constant. This is the simplest, non-pathological matter model with dynamical degrees of freedom in 
spherical symmetry. By this we mean the following: in spherical symmetry, Birkhoff's theorem completely 
rS ■ determines the local structure of electro-vacuum spacetimes, leaving no dynamical degrees of freedom; 

on the other hand, dust, for instance, is known to develop singularities even in the absence of gravity, 
i.e. in a fixed Minkowski background, and consequently is deemed pathological^ The self-gravitating 
scalar field appears then as an appropriate model to study gravitational collapse. This is in fact the 
original motivation behind the monumental body of work developed by Christodoulou concerning self- 
gravitating scalar fields with vanishing cosmological constant jj and it is inspired by these achievements 
that we proceed to the positive A casqj. 



1 For instance, either by symmetry conditions or smallncss assumptions on the initial data the formation of (cosmological) 
black holes is excluded in all the referred results. 

2 It should be noted that the presence of a positive cosmological constant may counteract the tendency of dust to form 
singularities. 

3 See the introduction to Chr09 for a thorough review of Christodoulou's results on spherically symmetric self-gravitating 
scalar fields. 

4 Christodoulou's work has also inspired a considerable amount of numerical work, including Choptuik's discovery of 
critical phenomena Cho93^ (see also IGM07| and references therein). The case A > seems to be less explored numerically, 
see however [Bra97 Bcy09b]. 
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We modify the framework developed in |Chr86j to accommodate the presence of a cosmological con- 
stant, thus reducing the full content of the Einstein-scalar field system to a single integro-differential 
evolution equation. It is then natural, given both the structure of the equation and the domain of the 
Bondi coordinate system where the reduction is carried out, to consider a characteristic initial value 
problem by taking initial data on a truncated null cone. 

For such an initial value problem we prove well posedness, global existence and exponential decay 
in (Bondi) time, for small data. From this, it follows that initial data close enough to de Sitter data 
evolves, according to the system under consideration, to a causally geodesically complete spacetime (with 
boundary), which approaches a region of de Sitter asymptotically at an exponential rate; this is a non- 
linear stability result for de Sitter within the class under consideration and can be seen as a realization 
of the cosmic no-hair conjecturc(f|. Also, we note that the exponential decay rate obtained, < e~ Hu , with 
H = 2y/A/3, is expected to be sharpd |Ren04j . Moreover, an interesting side effect of the proof of our 
main results is the generalization, to this non-linear setting, of boundedness of the supremum norm of 
the scalar field in terms of its initial characteristic data. We refer to Theorem [T] for a compilation of the 
main results of this paper. 

As was already clear from the study of the uncoupled case [CAN 12] . the presence of a positive 
cosmological constant increases the difficulty of the problem at hand considerably. In fact, a global 
solution for the zero cosmological constant case was obtained in |Chr86j by constructing a sequence 
of functions which, for an appropriate choice of Banach space, was a contraction in the full domain; 
such direct strategy does not work (at least for analogous choices of function spaces) when a positive 
cosmological constant is considered, since a global contraction is no longer available even in the uncoupled 
case. Moreover, new difficulties appear in the non-linear problem when passing from zero to a positive 
cosmological constant: first of all, the incoming light rays (characteristics), whose behavior obviously 
depends of the unknown, bifurcates into three distinct families, with different, sometimes divergent, 
asymptoticsQ this is in contrast with the A = case, where all the characteristics approach the center of 
symmetry at a similar rate. Also, for a vanishing cosmological constant the coefficient of the integral term 
of the equation decays radially, which is of crucial importance in solving the problem; on the contrary, 
for A > such coefficient grows linearly with the radial coordinate. 

To overcome these difficulties we were forced to differ from Christodoulou's original strategy consid- 
erably. The cornerstone of our analysis is a remarkable a priori estimate, the aforementioned result of 
boundedness in terms of initial data, whose inspiration comes from the uncoupled case [CAN12] . We can 
then establish a local existence result with estimates for the solution and its radial derivative solely in 
terms of initial data and constants not depending on the time of existence, which allows us to extend a 
given local solution indefinitely. The decay results, which in the vanishing cosmological setting are an 
immediate consequence of the choice of function spaces and the existence of the already mentioned global 
contraction, here follow by establishing "energy inequalities", where the "energy function" is given by the 
supremum norm of the radial derivative of the unknown (|66p . 

To make this strategy work we were forced to restrict our analysis to a finite range of the radial 
coordinate; one should note nonetheless, that although finite, the results here hold for arbitrarily large 
radial domains. At a first glance one would expect the need to impose boundary conditions at r = R, 
for R the maximal radius; this turns out to be unnecessary, since for sufficiently large radius the radial 
coordinate of the characteristics becomes an increasing function of time, and consequently the data at the 
boundary r — R is completely determined by the initial data (see Figure [2]). This situations parallels that 
of |Rin08) , where local information in space (here, in a light cone) allows to obtain global information in 
time. 

A natural consequence of the introduction of a positive cosmological constant is the appearance of 
a cosmological horizon. In fact, although the small data assumptions do not allow the formation of a 
black hole event horizon, a cosmological apparent horizon is present from the start, and a cosmological 
horizon formed; this is of course related to the difficulties mentioned above concerning the dynamics of 

5 Albeit in a limited sense, since our coordinates do not reach the whole of future infinity (see Figure [T{. A precise 
statement of the cosmic no-hair conjecture can be found in Bcy09a , where it is shown that it follows from the existence of 
a smooth conformal future boundary 

6 Although our retarded time coordinate u in Q is different from the standard time coordinate t, it coincides with t along 
the center r = 0, and hence is close to t in our r-bounded domain, thus giving the same exponential decay. For instance, in 

de Sitter spacetime we have u = t — y3/A log ( 1 + y/A/3 r j . 



7 The use of double null coordinates (u, v), also introduced by Christodoulou for the study of the Einstein-scalar field 
equations in Chr9F, would facilitate the handling of the characteristics, which in such coordinates take the form v = const., 
but in doing so we are no longer able to reduce the full system to a single scalar equation. 
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the characteristics. 



1.1 Previous results 

A discussion of related results in the literature is in order. 

The first non-linear stability result for the Einstein equations, without symmetry assumptions, was the 
non-linear stability of de Sitter spacetime, within the class of solution of the vacuum Einstein equations 
with positive cosmological constant, obtained in the celebrated work of Friedrich [FH86 1 This result 
is based on the conformal method, developed by Friedrich, which avoids the difficulties of establishing 
global existence of solutions to a system of non-linear hyperbolic differential equations, but seems to 
be difficult to generalize to Einstein-matter systems. A new, more flexible and PDE oriented approach 
was recently developed by Ringstrom |Rin08) to obtain exponential decay for non-linear perturbations of 
locally de Sitter cosmological models in the context of the Einstein-nonlinear scalar field system with a 
positive potential; these far-reaching results have a wide variety of cosmological applications, but break 
down exactly in the situation covered here, since they assume that the potential V satisfies V"(0) > 
(and so cannot be constant). 

In the meantime, based on Ringstrom's breakthrough, Rodnianski and Speck IRS091 . and later 
Speck |Spell| , proved non-linear stability of FLRW solution with flat toroidal space within the Einstein- 
Euler system satisfying the equation of state p = c s p, < c s < 1/3; exponential decay of solutions close 
to the flat FLRW was also established therein. In the same context, by generalizing Friedrich's conformal 
method to pure radiation matter models, Liibbe and Kroon |LK11] were able to extend Rodnianski and 
Speck's non-linear stability result to the pure radiation fluids case, c s = 1/3. 



1.2 Main results 

Our main results may be summarized in the following: 

Theorem 1. Let A > and R > \/3/A. There exists eo > 0, depending on A and R, such that for 
0o € C fc+1 ([0,i?]) (fc > 1) satisfying 

sup \<f>o(r)\ + sup \d r <f> Q (r)\ < e , 

0<r<R 0<r<R 

there exists a unique Bondi-spherically symmetric C k solution^ (M, g, <ft) of the Einstein- A-scalar field 
system Q, with the scalar field 4> satisfying the characteristic condition 

0U=o = <£o ■ 

The Bondi coordinates for M have range [0, +00) x [0, R] x S 2 , and the metric takes the form (TT|). Moreover, 
we have the following bound in terms of initial data: 

\(j)\ < sup \d r (rcj) (r))\ . 

0<r<R 

Regarding the asymptotics, there exists € R such that 

W,r)-±\ <e~ Hu , 

and 

where H := 2y / A/3 and g is de Sitter's metric in Bondi coordinates, as given in @. Finally, the spacetime 
(M, g) is causally geodesically complete towards the future^ and has vanishing final Bondi mas 

8 This was later generalized to n + 1 dimensions, n odd, by Anderson | And05| . 

9 See Section [2] for the precise meaning of a C 1 solution of the Einstein- A-scalar field system in Bondi-spherical symmetry. 
10 A manifold with boundary is geodesically complete towards the future if the only geodesies which cannot be continued 

for all values of the affine parameter are those with endpoints on the boundary. 
11 See Section [3] for the definition of the final Bondi mass in this context. 
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This result is an immediate consequence of Proposition [T] and Theorems |3] and 0] Note that, as is the 
case with the characteristic initial value problem for the wave equation, only <p needs to be specified on 
the initial characteristic hypersurfac^H (as opposed to, say, <fi an d d u 4>). There is no initial data for the 
metric functions, whose initial data is fixed by the choice of <f>Q. A related issue that may cause confusion 
is that the vanishing of <9 r 0o(O) is not required to ensure regularity at the center: in fact, the precise 
condition for to be regular at the center is d u <fi(u, 0) = d r <p(u, 0), which is an automatic consequence of 
the wave equation ([7])- The reader unfamiliar with these facts should note that, for example, the unique 
solution of the spherically symmetric wave equation in Minkowski spacetime, <9 2 (r</>) — d 2 (r<fi) = 0, with 
initial data 4>{r, r) — r, is the smooth function cf>(t, r) = t for t > r. 



2 Einstein-A-scalar field system in Bondi coordinates 

We will say that a spacetime (M, g) is Bondi- spherically symmetric if it admits a global representation 
for the metric of the form 

g = —g(u, r)g(u, r)du 2 — 2g(u, r)dudr + r 2 dQ 2 , (1) 

where 

dn 2 = d6 2 + sin 2 6dip 2 , 
is the round metric of the two-sphere, and 

(u, r) £ [0, U) x [0, R) , U, R £ M+ U {+00} . 

If U or R are finite these intervals can also be closed, which corresponds to adding a final light cone 
{u = U} or a cylinder [0, U) x S 2 as a boundary, in addition to the initial light cone {u — 0}; the metric 
is assumed to be regular at the center {r = 0}, which is not a boundary. 

The coordinates (it, r, 9, tp) will be called Bondi coordinates. For instance, the causal future of any 
point in de Sitter spacetime may be covered by Bondi coordinates with the metric given by 

g = - (l - y 2 ^j du 2 - 2dudr + r 2 dQ 2 (2) 

(see Figure [I}. Note that this coordinate system does not cover the full dc Sitter manifold (which 
strictly speaking is not Bondi-spherically symmetric), unlike in the asymptotically flat A = case. The 
boundary of the region covered by Bondi coordinates, that is, the surface u = — 00 (which for A = 
would correspond to past null infinity), is an embedded null hypersurface (the cosmological horizon of 
the observer antipodal to the one at r — 0). Moreover, the lines of constant r, which for A = approach 
timelike geodesies asr-> +00, here become spacelike for sufficiently large r. 

Although the causal structures of Minkowski and de Sitter spacetimes are quite different, the ex- 
istence of Bondi coordinates depends solely on certain common symmetries. More precisely, a global 
representation for the metric of the form ([T]) can be derived from the following geometrical hypotheses: 

(i) the spacetime admits a 5*0(3) action by isometries, whose orbits are either fixed points or 2-spheres; 

(ii) the orbit space Q = M/SO(3) is a 2-dimensional Lorentzian manifold with boundary, corresponding 
to the sets of fixed and boundary points in M; 

(iii) the set of fixed points is a timelike curve (necessarily a geodesic), and any point in M is on the 
future null cone of some fixed point; 

(iv) the radius function, defined by r(p) := Area,(O p ) / 4n (where O p is the orbit through p), is mono- 
tonically increasing along the generators of these future null cones 

12 Notice however that uniqueness is not expected to hold towards the past. 

13 These two last assumptions exclude the Nariai solution, for instance, from our analysis. 
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u = — oo 



Figure 1: Penrose diagram of de Sitter spacetime. The dashed lines u = constant are the future null 

cones of points at r — 0. The cosmological horizon H corresponds to r = and future infinity to 
r = +00. 

2.1 The field equations 

The Einstein field equations with a cosmological constant A are 

Rpv - -^Rg^is + ^9nv = kT^ , (3) 

where is the Ricci curvature of g, R the associated scalar curvature and T^ v the energy- momentum 
tensor. For a (massless) scalar field <f> the energy-momentum tensor is given by 

= d^dvt - \ g ap (d a <t>) (dp<i>) , 

and then ([3]) becomes 

R» v = nd^(j)d u (j) + Ag^ v . (4) 

These equations are written for a spacetime metric of the form (flj and a spherically symmetric scalar 
field in the Appendix. As shown in [Chr86] . their full content is encoded in the following three equations: 
the rr component of the field equations, 

"a = \OrY) ; (5) 

r g or 

the 89 component of the field equations, 

|-(rg)=.g(f-Ar 2 ) ; (6) 
and the wave equation for the scalar field, 

v% = & va^ = , (7) 

which reads 

lid g d 

r du 2 dr dr v ' r ' 2 \ dr ) \ dr 



8 w>-Ug'ii£i- w 
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2.2 Christodoulou's framework for spherical waves 

Integrating ([5]) with initial condition 

g{u,r = 0) = 1 

(so that we label the future null cones by the proper time of the free-falling observer at the centeiF^) 
yields 

g = e t i;<a^) 2 ds_ (g) 

Given any continuous function / = /(m, r) we define its average function by 

f(u,r):=- ( f(u,s)ds, (10) 
r Jo 

for which the following identity holds: 

Or r 

Using the regularity condition 

lim rq = , 

implicit in our definition of Bondi-spherically symmetric spacetime, we obtain by integrating (pj): 

1 r r A r r 

g = - g(l- As 2 ) ds =g(l- Ar 2 ) = g / gs 2 ds . (12) 

r Jo r Jo 

Following |Chr86j we introduce 

h := d r {r<t>) . 

Assuming (f> continuous, which implies 

lim r<f> = , 

r— ¥0 

we have 

, 1 r , , s , r ,dct> dh h-h 

<b — — h (u, s) ds = h and — = — = , (13) 

r J Q or or r 

and so ^ reads 

5 ( M ,r) = cxp (j^ dsj . (14) 

Now, defining the differential operator 

D ■= — - I — 

du 2 dr ' 

whose integral lines are the incoming light rays (with respect to the observer at the center r = 0), and 
using (|13[) together with ([B]), the wave-equation ((5J) is rewritten as the integro-differential equation 

Dh = G(h-h) , (15) 

where we have set 

G := l -d r g (16) 
1 

27 



(1 - Ar 2 )g - (1 - Ar 2 )g (17) 



(g-g) A f r 2 A 

— / gs ds - -rg . (18) 



2r 2r 2 J " 2 

Thus we have derived the following: 

Proposition 1. For Bondi-spherically symmetric spacetimes ([TJ, the Einstein-scalar field system with 
cosmological constant (HJ) is equivalent to the integro-differential equation (fT5)l . together with (TT2|) . (fT5|) . 
(HU) and (HSJ. 



14 This differs from Christodoulou's original choice, which was to use the proper time of observers at infinity. 
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We will also need an evolution equation for d r h given a sufficiently regular solution of (|15l) : using 



[D,d r ]=Gd r , 

differentiating (1151) . and assuming that we are allowed to commute partial derivatives, we obtain 

Dd r h-2Gd r h = -Jd r h, (19) 

where 

J: = G-rd r G (20) 



= 3G + Agr + (Ar 2 - 1) 



ldg 
2dr 



(21) 



3 The mass equation 

Consider a Bondi-spherically symmetric C k solution of Q on a domain (u, r) £ [0, U) x [0, R] (with 
R > y/3/A). From equations (JSJ) and (|14p it is clear that rg is increasing in r for r < y/l/A and 
decreasing for r > v/l/A. On the other hand, equation (|12[) implies that <jp(u, r) approaches — oo as 
7' — > +oo. Therefore there exists a unique r = r c (u) > y/l/A where g(u,r) vanishes. This defines 
precisely the set of points where -j^ is null, and hence the curve r — r c (u) determines an apparent 
(cosmological) horizon. Since g is increasing in r, we have from (|12p 



, r) < g(u, VW)\f Q (1 - As 2 ) ds = g(u, jl/A) (l - ^ 



Therefore the radius of the apparent cosmological horizon is bounded by 

x < rc{u) - n 

for all u. From ([6]) it is then clear that || < for r = r c (u), and so by the implicit function theorem the 
function r c (u) is C k . From the uu component of H]) (equation ([75)1 in the Appendix), we obtain 



r 9m \.g 



n(d u (t>f 



when g = 0, showing that ^ is nondecreasing in u, and so r c (u) must also be nondecreasing. In particular 
the limit 

n := lim r c (u) 

U-¥U 

exists, and v/l/A < ri < y/Z/A. 

We introduce the renormalized Hawking mass functional |Nak95l IMN08j 



m(w, r) 



A 
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(22) 



which measures the mass contained within the sphere of radius r at retarded time u, renormalized so as 
to remove the contribution of the cosmological constant and make it coincide with the mass parameter 
in the case of the Schwarzschild-de Sitter spacetime. This function is zero at r = 0, and from (O, © we 
obtain 

dm nr 2 g 2 

implying that m(u, r) > for r < r c {u). We have 

r c (u) 



m(u,r c (u)) 



1 A /- \2 



15 This function is also known as the "generalized Misner-Sharp mass". 
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whence 

±m(u, r c {u)) = tM- (l - Ar c (u) 2 ) < 0, 
and so m(it, f c (w)) is a nonincreasing function of u. Therefore the limit 

Mi := lim m(u,r c (u)) = ^ ( 1 - ^-r 2 ) 

exists, and from v/l/A < rx < -\/3/A we have < Mi < 1/V9A- We call this limit the /maZ Bondi mass. 
Note that, unlike the usual definition in the asymptotically flat case, where the limit is taken at r — +oo, 
here we take the limit along the apparent cosmological horizon; the reason for doing this is that r < R 
in our case. 



4 Basic Estimates 

Given U,R>0, let Cfj R denote the Banach space (c°([0, U] x [0, R]),\\- || c ° H ) , where 



l/ll C o R := 

U ' R (u,r)e[0,C7]x[0,.R] 



sup |/(w,r)|, 



and let denote the Banach space of functions which are continuous and have continuous partial 

derivative with respect to r, normed by 

\\f\\xo, R ■= \\f\\c° iR + \\ d rf\\c« UR ■ 

For functions defined on [0, R] we will denote C°([0,R]) by C R , C 1 ([0, R}) by Xr, and will also use these 
notations for the corresponding norms. 



For h € C[Jr R we have 



1 



h(u,r)\<- \h(u,s)\ds < - \\h\\ c o ds = \\h\\ c 



and if /i € Xr/ b we can estimate 



(/i — h)(u,r)\ 



(h(u, r) — h(u, s)) ds 


r />r 



dh 
dp" 



(u, p)dp ds 



< 



1 r r r 

r Jo J " 9r/l " Cg -.« dpdS = 2 " 9r/l " c &,fl 



(23) 



Thus 



and by (fT4")l we get 



2 7 



dr< ^\\d r hf c0 R\ 
16 >. 



,/( „. 0) = 1 < .9(71, r) < X := exp ( ^R/i||*o „i? 2 



(24) 



4.1 The characteristics of the problem 



The integral curves of D, which are the incoming light rays, are the characteristics of the problem. These 
satisfy the ordinary differential equation, 



dr 
du 



1 



g(u,r). 



(25) 



To simplify the notation we shall denote simply by u H- r(u), the solution to (l25l) . satisfying r(u\) = r\. 
However it should be always kept in mind that r(u) = r(u;ui,ri). 
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Using (J23| we can estimate g, given by (fT2")l . and consequently the solutions to the characteristic 
equation (j2"ij)) : In fact, for r < -j- => 1 — Ar 2 > we get 

§>- T(l - As 2 )ds = 1 - -r 2 > 1 - K-r 2 . 
r Jo 3 3 

For r > -7= we have 

— VA 



g(u, r) > i /"^ (l - As 2 ) ds + ^- J (l - As 2 ) ds 



2 (l-X) + xfl--r 2 



>|(l_f0+*(l-§r 2 
We then see that the following estimate holds for all r > 0: 



3 

The same kind of reasoning also provides the upper bound 



KA n 

9 > 1 - ^rr 2 . (26) 



3 

From ([23)1 and (|2"6"|) we now obtain the following differential inequality 

dr 1 AJf 



g<K--r 2 . (27) 



Denoting 



a = —\ and r„ = \ - — 

2 V 3 V AK 



where r c is the positive root of the polynomial in (|26[) . the solution r (u) of the differential equation 
obtained from (|28p (by replacing the inequality with an equality) satisfying r~ (ux) = r\ < r~ is given by 

r~ (u) = — tanh ja(c~ — it)} , 

for some c~~ = c _ (ui, n); by a basic comparison principle it then follows that whenever r(u±) = r\ < r~ 
we have 

r(u) > — tanh |a(c~ - u)\ , Vu < u x . (29) 
2a 

Denote the positive root of the polynomial in (|2"T|) by 

then, for appropriate choices (differing in each case) of c~ = c _ (ui,ri) and c + = c + (ui,ri), similar 
reasonings based on comparison principles give the following global estimates for the characteristics (see 
also Figure HJ) : 

• Local region {r\ <r~): 

— tanh |a(c~ — u)\ < r(u) < — tanh |a(c + — u)j , Vu < u\ . (30) 
2a 2a 

• Intermediate region (r~ <r\ < r~£~): 

1 K 

— coth {a(c~ — u)} < r(u) <— tanh |a(c + — u)\ , Vw < u\ . (31) 
2a 2a 

• Cosmological region (r > r+): 

1 if 

— coth {a(c~ - u)} < r(u) < —coth {a(c + - u)} , Vu < u x . (32) 
In particular, for r(iti) = r\ > r~ we obtain 

r(ti) > r" > , Vu < mi. (33) 
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Figure 2: Bounds for the characteristics through the point (ui,r\) in the local (r\ < r c ), intermediate 
{ r c — r i < r c") an d cosmological (ri > r+) regions. 

4.2 Lemma 1 

The purpose of this section is to prove the following lemma: 

Lemma 1. Let A > and R > 0. There exists x* = x*(A,R) > and constants C% — Ci(x* , A, R) > 0, 
such that if \\h\\ Xu R < x* , theiQ 

G<~C ir , C\ = ^ + 0(x*), (34) 

\G\ < C 2 r , (35) 
|J|<C 3 r, C 3 = 0(x*), (36) 

and, for any u\ > and r\ < i?, 

"1 / fll \ 

exp f y 2G(u, r(v))dvj du < C 4 , (37) 
where r(u) = r(u; u\, r\) is the characteristic through (u\, r\). 

Remark 1. We stress the fact that while allowed to depend on R the constants do not depend on any 
parameter associated with the u-coordinate. 

Proof. We have, from (J23J), 

1<9<K* :=cxp^(x*) 2 R 2 ) • (38) 
Differentiating (fJU) while using ([23]) and (J3SJ leads td3 

< d JL < gil^R- < K * (x*)\ (39) 



16 As usual, 0(x*) means a bounded function of x* times x* in some neighborhood of x* = 0. 

17 From now on we will use the notation / < g meaning that / < Cg, for C > only allowed to depend on the fixed 
parameters A and _R. 
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and consequently 



< (g - g) (u, r) = - I (g(u, r) - g(u, s)) ds 



< 



-II K*(x*) 2 pdpds 

r JO Js 



< K*{x*) 



From this estimate, (|18p and (|38p we see that 
'A A 



K* r < G < 



^* c{x*y + 



for some constant C > depending only on A and i?. Since K* — > 1 as x* 
choosing x* appropriately small. Also, inequality (|35p is immediate. 
From ([21]), QS) and (gDJ) we now obtain ([36) . 



To prove (|37|) we start by using ([34)) to obtain 



j/" 1 2G(v,r(v))dv du 



< 



If n < r c = W then ([2"§)> holds and we then have 



Since 



and 



we obtain 



rill q rui 

-2Ci / r(v)dv < / tanh (a(cT - v))dv 

Ju a Ju 

C\ j / cosh (a(c~ — Mi)) 
a 2 \ cosh(a(c~ — it)) 

cosh(a(c- -m)) < 2eQ;( „_ Ul ) 
cosh (a(c _ — u)) — 



1 /A 1 /AX . /— 

-\ — < a = -\ < VK* 

2V 3 " 2 V 3 ~ 



(40) 

0, ([Ml) then follows by 



Ci 

as desired. 

If ri > , we have ([33]) which gives 



2 Cl J^r(v)dv du< / e -2C 1 r e -(u 1 -«) du< 



L "1 



< 2 Ci/q2 -^ < C 4 (a;*,A,i?) 



(41) 



2Cir c - 



1 - e 



-2Cir~m 



2(T 

<7f <C4(x*,A,i?) 



which completes the proof of the lemma. 



□ 



5 Controlled local existence 

Local existence will be proven by constructing a contracting sequence of solutions to related linear prob- 
lems. Given a sequence {h n } we will write g n := g(h n ), G n :— G(h n ), etc, for the quantities ([T4]) . ([To]) , 
etc, obtained from h n ; for a given h n the corresponding differential operator will be denoted by 

D n = d u - , 
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and the associated characteristic through (wi,ri) by Xn — Xn(u) — (u, r n (u; Ui, as before, we will 
drop the explicit dependence on initial conditions when confusion is unlikely to arise. 

With these notational issues settled we are ready to prove the following fundamental result: 



Lemma 2. Let A > 0, R > Jf and h G ^([O,^]). There exists x* = x*(A,R) > and C* 
C*(x*,A,R) > such that if 

x* 

\\ho\\x R < Y^~c* ' 
then the sequence {/i„}„gN„ defined by ho(u, r) = hg{r) and 

D n h n +i — G n h n +\ = —G n h n 
h n+ i(0,r) = h (r) , 



is in C\[0, R] x [0, U]) and satisfies 



G n < 0, (42) 
IIMc° = IIMc°, (43) 
< (l + C*)\\h \\ XR , (44) 



■'U.R 



for all neNo and all U > . 

Remark 2. We stress the fact that C* does not depend on either U or n. 

Proof. The proof is by induction. That the conclusions follow for the th term is immediate, with (l42j) 
obtained from Lemma [T] by setting x* accordingly small. Assume that h n satisfies all the conclusions of 
the lemma. In particular, since we have h n e C ll ([0, U] x [0, R}) we see, from the respective definitions, 
that h n , g n and g n are C 1 for r ^ 0; regularity at the origin then follows by inserting the first order 
Taylor expansion in r of h ni centered at r = 0, in the definitions of h n , then g n and finally g n . Later in 
the proof we will also need d r G n to be well defined and continuous in the domain under consideration; 
this follows by using the previous referred expansions in equations (|20[) and (|21[) . 

Note that, as a consequence of the regularity for g ni we also obtain well posedness and differentiability 
with respect to the initial datum r\ for the characteristics given by (|25p ; in particular we are allowed to 
integrate the linear equation for h n +i along such characteristics to obtain 

K+xiuun) = /icMO))^' 1 G ^ dv - f 1 (G n hn). e^ 1G ^ dv du . (45) 



This defines a function h n+ i : lZ n +i C [0, U] x [0, R] — > K where 

Ttn+i = {(u,r) | Xn(u) = (u,r n (u)) = (u,r) and r„(0) G [0,R]} . 

Since the problem for the characteristics is well posed, there is a characteristic through every (u, r) G 
[0, U] x [0, R]; in particular lZ n +i is non empty, but nonetheless, integrating backwards in it, the char- 
acteristics may leave the fixed rectangle before reaching u = 0, which in turn would lead to 1Z n +i 7^ 
[0, U] x [0, R]. We may rule out this undesirable possibility by a choice of appropriately small x*; in fact, 
it suffices to guarantee that the r„ component of all characteristics with sufficiently large initial datum r% 

are nondecreasing in u: given R > v/^, since (|44[) and the smallness condition on the initial data imply 

\\hn\\xu, B < (l + C*)\\h \\x R <X* , 

we see that (recall (|24|l ) 

K n < K* - e c ^ R2 , 
and from the global characterization (|5U|) - (l52l the desired monotonicity property follows if 



K n A 

-l-<R, 

which can be arranged by choosing x* sufficiently small (see also Figure [2]). 
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We have already showed that h n , Xn and G n have continuous partial derivatives with respect to r; 
we are then allowed to differentiate (I45j) with respect to r\ and, since D n h n+ \ is clearly continuous, we 
conclude that 

h n+ i E C^iO,!/} x [0,R]) . 

From the previous discussion d r Dh n+ i is continuous, so differentiating equation (|4"5j) with respect to r, 
and using the fact thatf^l 

[D n ,d r ] = G n d r 

we obtain the following differential equation for d r h n+ i (recall ([20)) ): 

D n (d r h n +i) — 2G n d r h n+ i = d r G n (h n +i — h n ) — G n d r h n 

_ _ dhn _ , t _ r , (K+i - K) 
or r 

Using the initial conditions 

d r h n+1 (0,r) = d r ho(r) 
and integrating along the characteristics leads to 



d r h n+1 ( Ul ,ri) = d r h (xn(0))e f ° l2G "^ dv 



[ i ah 4- r T -G \ ~ h n) 



lx» 



By the induction hypothesis we have G n < and ||ft>n||c — ll^ralle — ll^ollc i therefore 



Then 



(46) 



|fcn+i(tti,ri)| < \\h \\ c o e J« lG ^ dv + fh n \\ c » / -G n \ Xn eJ" 1G ^ dv du 

' Jo 

< \\ho\\c° R (e f " G ^ dv - jT G n \ Xn e^ G ^« d *dv\ = \\h \\ e o R 



\h — h I 2||/in||r-o 

\h n+1 -h n \<2\\h Q \\c% and \d r h n \= 1 -^ ^ < (47) 



so that, relying once more on Lemma [TJ 

\d r h n+1 ( Ul , ri )\ < \\d r ho\\ c0 efo 1 ^\ Xn dv 

+ 2{C 2 + 2<7 3 )IMIcS r e/:i 2GlAxndvdu 
R Jo 

< \\d r h \\ c0 +2(C 2 +2C 3 )C 4 \\h Q \\ c a . 
Setting C* := 2(C 2 + 2C 3 )C 4 it now follows that 

\\hn+i\\xu, R = \\h n +i\\c° R + \\d r h n +i\\c° iR 

< (l + C*)||/io||c° + RMc° 

< (l + C*)\\h \\ XB . 

Thus, if || h \\ Xr < if^r, then 

\\ h n+l\\x U:R < X* 

and by Lemma [1] 

G„+i < , 

which completes the proof. □ 



18 Here we are using the following generalized version of the Schwarz Lemma: if X and Y are two nonvanishing C 1 vector 
fields in R 2 and / is a C 1 function such that X ■ (Y ■ f) exists and is continuous then Y ■ (X ■ /) also exists and is equal to 
X-(Y-f)-[X,Y]-f. 
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Lemma will now allow us to establish a local existence theorem for small data, while controlling the 
previously defined supremum norms of the solutions in terms of initial data. 

Theorem 2. Let A > 0, R > and h £ C k ([0,R}) for k > 1. There exists x* = x*(A,R) > and 
C* (x* , A, R) > 0, such that, if ||/io||x n < TTc* » then the initial value problem 



Dh = G(h-h) 
h(0,r)=h (r) 



(48) 



has a unique solution h G C fc ([0, C/] x [0, i?]), for U = U(x*/(1 + C*)\ R, A) sufficiently small. Moreover, 



and 



\h\\ XuR <(l + C*) \\h \\ XR 



(49) 
(50) 



Proof. Fix x* as in Lemma[T]and consider a sequence {h n } as defined in Lemma[3J with U < 1. From 
and Lemma [5] we have 

\(h n - h n ) + {h n -x - h n -i)\ < T -z (||#rMlc° r + \\d r h n -i\\ c o r J 

< (l + C*)r\\h \\x R <x*r, 

and 



|(/i„ - h n ) - {hn-x - /i„_i)| = (/i„ - /i„-i) - (h n - h n ~i) < 2 \\h n - 7i n _i|| c o 

so that 

\(h n - h„) 2 - (/i„_i - h n ^i) 2 \ < \ (h n - h n ) + (h n -i - h n -i)\ \{K - h n ) - (/i n _i - h n -i)\ 

< 2x* r \\h n - h n -i\\ c o 

= Cr \\h n - h n -i\L 



(51) 



(we will, until the end of this proof, allow the constants to depend on x* , besides the fixed parameters A 
and R) . The mean value theorem yields the following elementary inequality 



\e x -e v \ <max{e x ,e^}|a;-j/| 



(52) 



from which (recall (|55|) ) 

\g n - #71-1 1 



exp C 



r {h n ~ h n f 



exp C 



{h n -i — h n -i) z 



< K* 



\(h n — h n ) 2 — (/ln-1 ~ hn-x) 2 



-ds 



(53) 



<Cr\\h n — h n -i\\c° 



Then 



and using (fTTJ) , 



in - yn-1 



(9n ~ 9n-i)(l - As 2 )ds 



\Gn ~ Gn-l\ — 7T~ 
Zr 



< Cr \\h n - h n -i\\ e o 



(ffn - Sn-i)(l - Ar 2 ) - - I (g n - g n -i)(l ~ As 2 )ds 



(54) 



< C \\h n — h n -\ c o 

Note that, since r < R, the r factors in the previous estimates may be absorbed by the corresponding 
constants. 
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Until now we have been estimating the difference between consecutive terms of sequences with both 
terms evaluated at the same point (u, r), but we will also need to estimate differences between consecutive 
terms evaluated at the corresponding characteristics; more precisely, for a given sequence f n we will 
estimate 

l/n|x« -/n-l|x»-il = \.fn(u,r n (u)) - /„_i(u,r„_i(u))| 

< \fn(u,r n (u)) - /„(u,r n _i(u))| + \f n (u,r„-i(u)) - f n -i(u,r n -i(u))\ . 

If for the second term we have, as before, a uniform estimate of the form C\\h n — /i n _i||c° , and for the 
first one of the form C\r n — r„_i|, then, by (|5rj|) below, we will obtain, since u\ < U < 1, 



\fn\xn ~ /n-l|x„_il - C\\f n - fn-l\\c^ >R 

Also, if ||9 r /„|| c o < C then the desired 



(55) 



\fn(u,r 2 ) - f n (u,ri)\ < / d r f n (r)dr < C\r 2 - r x \ , 

Jri 

follows immediately. We have (see 

\d r h n \= <C, 
r 

and from ([M| 

\d r g n \ < Cr . 

By Lemma[T]we have ||Gri||e° R < C, which in view of (TTO)) is equivalent to ||3r<7n|leg. — @\ since (f^U|) . d2U) 
and ([55)1 together with the above bounds yield ||<9 r G n || c o ^ < C, the desired estimates, of the form (jS"5")) . 

follow for the sequences h n , h n , g n , g n and G n once we have proved (|56p . To do this, start from 
equation ((25|) for the characteristics to obtain 



r n (u) = r n (ui) + - J g n (s,r n (s))ds , 

so that the difference between two consecutive characteristics through (ui,ri) satisfies 
1 f Ul 

r n {u) - r n -i(u) = - / {<7„(s, r„(s)) - g n _i(s, r n _i(s))} ds 



{<7„(s,r„(s)) - 9n (s, r„_i (s) ) } rfs + - / {g n (s, r„_i (s)) - <?„-i(s, r„_i (s))} ds . 



From the previously obtained bounds ||9 r ffn||c° — C an d (|54|) . we then have 



\r n {u) - r n _i(u)| < C / \r n (s) - r n -i(s)\ds + C (ui - u)\\h n - h n -i\\ c o 



from whicl 
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|r n («) - r„_i(u)| < - ^-i|lc«, fl (e c <— > - l) 



(56) 



as desired. 

Now, from f)45[) and the elementary identity 



a 2 b 2 C2 - a\b\C\ = (a 2 - ai)6 2 c 2 + (62 - &i) a i c 2 + ( c 2 - ci)ai^i 



19 Here we used the following comparison principle: if y,z £ C°([to,*l]) satisfy j/(t) < f(t) + C Jl 1 y(s)ds and 2(i) 
/CO + C ft 1 z(s)ds, then j/(t) < z(t), Vt G [to,*l] ■ 
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we get 

\(h n+1 - h n )(ui,n)\ < \\h \\ c o 



exp ( / G n \ Xn dv ) — exp ( / G 



-i|x»- 



dv 



L 



+ / | G «|X~ _G «-l|x™-i| r«lxJ eX P / G n\ X n dv ) du 



II 



\h n \Xn - K-l\ Xn -l | \ G n-l\xn-l I eX P / G n\ Xn dv du 



III 



exp ( / G n i x dv I — exp ( / G 



T n-l\Xn- 



, dv 



l^n-llXn-i^n-llXn- 



Using dHJ), (15"2"j) and (|55|) . which holds for the sequence G n as discussed earlier, gives 



G n \ Xn dv 



G 



n-l|Xn- 



, dv 



< Cu\ \\h n - /l„_i|| c o 



7 < 

and, in view also of ((43]), 

II < C\\h \\c° \\h n - /in-lLo C?li < Cmi ||/l„ - /l„_i|Lo 

In a similar way (recall that (|55|) also holds for the sequence h n ) 

IIIKCutWK-hn-t 

and, using the bound for /, 



\c° i 



7F<Cu?||/i n -/i n _i|| c o . 
Putting all the pieces together yields (recall that we have imposed the restriction u\ < U < 1) 

\\h n+ i - h n \\ c a r < CU\\h n - fo n _i|| c o. R . 
Now, applying the same strategy to (pf6")l leads to 



(57) 



\(d r h n+1 - <9 r /i„)(wi,ri)| < ||d r /io||c° 



exp 2G„| Xri dv^ — exp 2G 



n-l|Xn- 



I J «lx„ ~ J n-i| x „_ 1 | \drK\ Xn \ exp / 2G n \ Xn dv \ du 



(ii) 

+ 1 \9 r h n \ Xn - d r hn-i\x n -i\ |^»-i|xn-i| ex P ( / 2G„| x „dw ) rfu 



(iii) 



rill 


exp ^ 


/■Mi 


Jo 







exp / 2G n \ Xn dv -exp / 2G„_i| Xn _ 1 cb 



r«-l|Xn-l^^n-l|Xn-l | dlL 



I 7 n G Ti 



|/Y 



n+1 'i-nl 



(*«) 

e f? 2G ntxn dv du 



\Xr, 



\Jn-l ~ G n-1 



(v) 

\h n — h n _i\ 
r 

(vi) 



IXn- 
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We have 

\d r g n - d r g n -i\ < 
< 



(h n — h n ) 2 (h n -i — h n _i) 2 

9n <7n-l 



9n 



r r 
\{h n — h n ) 2 — (h n -x — h n -i) 2 \ 



+ \dn ~ 9n-l\ 



(h n -i — h n -iY 



where we have used 



< \\hn - hn-l\\c° )R , 

(|5T|) and $55$. Similarly 



\d r g n (u,r 2 ) - d r g n {u,ri)\ < \g n (u,r 2 )\ 



(h n - h n ) 2 (u, r 2 ) (h n -h n ) 2 (u,ri) 



r 2 n 

\g n (u,r 2 ) -g n (u,ri)\ 



+ 



(/i„_i - /i„_i) 2 (m, ri) 



< 



\r 2 - ri\ 



and (|57p provides 
We finally obtain 



We conclude that ([55]) holds for the sequence d r g n and since it also holds for the sequences g n and G„ 
we obtain from ([2"T]) 

\Jn\Xn ~ •Ai-llxn-il - Cll^n _ ^n-lllcg, R ■ 

As an immediate consequence one obtains for (i) — (iv) estimates similar to the ones derived for / — IV 
(recall CD}, O and 

Using (|55|) and (155)) we also have 

< CU\\h n 
(v)<CU 2 \\h n -h n ^\\ c o R . 

\\h n +i — h n \\x a ,R = \\h n +i — h n \\ c o + \\d r h n+ i — d r h n \\ c o 

< CU\\h n - h n -x\\ c a R 

< CU\\h n - h n -i\\xu,R ■ 

So, for U sufficiently small, {h n } contracts, and consequently converges, with respect to || • Hx^r- The 
previous estimates show that the convergence of h n lead to the uniform convergence of all the sequences 
appearing in (05]) and (US]). Taking the limit of (fl5l) leads to 

h( Ul ,n) = fc o (x(0))e^' 1 G ^ d - - P (Gh) \ x ef? G ^ dv du , (58) 

Jo 

where we denote the limiting functions by removing the indices. Equation (|58[) shows that h is a continuous 
solution to (|48[) . the limit of (HHJ) shows that d r h solves (JT5J) and is continuous, and we see that h € C 1 , 
since Dh is also clearly continuous. 

Now let 1 < to < k be an integer, and assume that h E C m . As in the proof of Lemma (2] but using 
the Taylor expansion of order to, we can show that h, g (which controls the characteristics), G and d r G 
are also C m , from which it follows that d r (Gh) is C m . Taking the partial derivatives of ([58)1 as in [CAN 12] 
(using the assumed regularity of the initial data) we then see that actually h s C m+1 , and so h € C . 

To establish uniqueness consider two solutions of P5|) and derive the following evolution equation for 
their difference: 

Dx(h 2 - hi) - G x (h 2 ~hi) = \ (92 ~ gi)d r h 2 + (G 2 - G x ) (h 2 - h 2 ) - d (h 2 - h) . (59) 



Integrating it along the characteristics associated to hi yields 



|(/i 2 -/ii)(ui,r 1 )| < 



~ \g 2 - 5x| \0 r h 2 \ + \G 2 - d| \h 2 -h 2 \ + \Gi\ \h 2 - h x \ 



j: 1 G n^du 



Ixi 
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Setting 

S(u) = ||(>2-/ii)(tv)||c° . 
then, arguing as in the beginning of the proof of this theorem, we obtain, from the previous inequality, 

8{ui) < C / S(u)e^ 1Gl ^ dv du . 
Jo 

Applying Gronwall's inequality we conclude that 

S(u) < , 

and uniqueness follows. 

The estimates and (|5U)) are now an immediate consequence of Lemma [5J □ 

6 Global existence in time 

Theorem 3. Let A > 0, R > Wf and h £ C k ([0,R}) for k > 1. There exists x* = x*(A,R) > and 
C* (x* , A, R) > 0, such that, if H/iollx^ < JTfcFJ 1 ' then the initial value problem 

[ h{0,r) = h {r) 
has a unique solution h £ C k ([0, oo] x [0,i?]). Moreover, 

||^||c°([0,oo)x[0,ii]) = IIMc°([o,-R]) > ( 61 ) 

and 

IWIjf([0,oo)x[0,.R]) < (1 + C*)||/i ||x([0,fl]) • ( 62 ) 

Also, solutions depend continuously on initial data in the following precise sense: if h 1 and h 2 are two 
solutions with initial data /ij and h^, respectively, then 

Wh 1 - h 2 \\ e o uR <c(u,R,A)\\hl-h 2 \\c R , 

for all U > . 

Proof. From Theorem [5] there exists a unique h 1 £ C k ([0,Ui] x [0,R]) solving (1501) . with existence time 
Ui = U(x*/(1 + C*) 2 ) . Moreover 



\h l {U u -)\\ Xr < Wh^x^ < (1 + C*)\\h \\x R < T 



C* 



So Theorem H] provides a solution h 2 £ C fc ([0, J7 2 ] x [0,i?]) with initial data /i 2 (0,r) = h l (JJ\,r) and 
existence time U 2 = U(x*/(1 + C*)). Now, 

h : [0, U x + U 2 ] x [0, R] ->• K 

defined by 

• \ /. 2 (u,r) , u£[U 1 ,U 1 + U 2 ] . 

is the unique solution of our problem in C fc ([0, U\ + U 2 ] x [0, i?]). Since applies to both h 1 and h 2 we 
see that: 

ll' 1 llcg. liS = \\ h o\\c° R > 

so that 

ll^llc° 2 , H = 11^(^1, Ollco < HMco , 

and hence 

IWk i+U2(fl = |Mle° . (63) 
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Arguing as in the proof of Lemma [U we see that d r Dh is continuous and consequently d r h solves (|T9"f so 
that: u 

d r h(u 1 ,r 1 ) = d r h ( X (0))e^ l2G ^ dv - [ 1 (Jd r h). e^ 12G ^ dv du. (64) 

Jo lx 

Consequently, 

\d r h( Ul , ri )\ < \d r h {r )\e^ l2Gdv + r \J\\d r h\e^ 12Gdv du 

Jo 

< \\d r h \\ c0 +2C 3 c<4||M c o 

^R C H 

< \\d r h \\ c0 +C*\\h Q \\ c0 , 

R '-R 

where we have used an estimate analogous to (|47|) , the fact that Lemma Q] applies to h (with the same 
notation for the constants), and the fact that we may choose C* :— 2(26*2 + 63)64, which can be traced 
back to the proof of Lemma O 

Combining the last two estimates with the smallness condition on the initial data leads to: 

\\h\\x Ul+U2 . R < (1 + C*)\\h \\x R < ■ (65) 

So, by Theorem[21 we can extend the solution by the same amount U2 = U(x*/(1 + C*)) as before; the 
global (in time) existence then follows, with the bounds (|6ip and (p32"j) a consequence of and (l6"5l) . 

The continuous dependence statement follows by applying Gronwall's inequality to the integral in- 
equality obtained integrating equation (|59p and using the estimates derived in the beginning of the proof 
of Theorem [2 

□ 



7 Exponential decay 

Theorem 4. Let A > 0, R > ^ 

h G C fe ([0,oo] x [0,R]), of dSOl satisfies 



Theorem 4. Let A > 0, R > J ^ and set H = 2y j. Then, for ||/io||x R sufficiently small, the solution, 

ki 



sup \d r h(u,r)\ < Ce- Hu , 

0<r<i? 

and, consequently, there exits ftgl such that 

\h(u,r)-h\ < Ce~ Hu , 
with constants C and C depending on ||ft-o||xfj, R and A . 
Proof. Consider the solution provided by Theorem [3] Set 

£{u):=\\d r h{u,-)\\ c o R , (66) 

and 

E(uq) := sup £ (u) . 

U>Uq 

Arguing as in (|2"5)l we get 

\(h-h)(u,r)\< r -£(u). (67) 

Lemma [T] applies and note that, for a fixed xq > 0, the estimates ([55)1 and (|5rj|) are still valid, with 
x* replaced with E(uo), for the functions G and J restricted to \uq,oo) x [0, R]. Integrating (TTHJ) with 
initial data on u — uq gives, for u\ > uq (compare with (l64|) ) 



£ (ui) < £(u Q )e- 2Cl C + / £ (u) e - 2 ^ 



'«0 

by using (joTl) and (THil) ; once again we have used the notation for the constants set by Lemma [T] 
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Recall that r(u) = r(u; Ui,ri) and that if r± < r c = y jj< then, as in the calculations leading to (|4Tj) . 
we have 

£(ui) < £(« )2 Cl/tt! e^" 1 + 2 c ^ a2 - 1 C 3 R / £{u)e^^ u - u ^du , 

so that u 

e^ Ul £( Ul ) < 2 Cl / a2 £(u Q ) + 2 Cl ' a2 ' 1 C z R [ ' £(u)e^ u du . 



Applying Gronwall's Lemma to J-(ui) := e <x Ul £(ui) then gives 

e^ ui £ («i) < 2 Cl / Q2 f (tt ) exp ^/""^CaiZfa - u )) , 

so that finally 

£(ui) < 2 Cl / Q2 £(u )exp|^2 Cl / Q2 - 1 C 3 i?- ^ uij . 
For r\ > r~ we have (|33[) instead and a similar, although simpler, derivation yields 

£(ui) < £(u )exp\ ( — |^ - 2Cir~ ) in 



Observe that if = e °( E <>°», Cl = f + 0(E(uq)), C 3 = O(E(u )), uniformly in u since u ^ £(uo) is 
bounded. Using such boundedness once more, we can encode the previous estimates into 

£{u) < Ce~ 6 ^ u , (68) 

with 

H(u ) = H + O(E(u )) . (69) 
Since E(uq) is controlled by ||/io||x B (see (J62J) ) , choosing the later sufficiently small leads to 

H(u ) >H>0, 

£{u) < Ce-" u 
E(u ) < Ce-" U0 , 
\H-H(u )\ < Ce- kuo . 



so that ([681) implies 
for u > uq. Then clearly 
so that (f6"9"|) becomes 
Finally, setting uq — % yields 



e Hu £{u) < Cexp(Hu - H(u/2)u) 
< Cexp(Ce- 6u ^u) < C , 

as desired; the remaining claims follow as in |CAN12j . □ 
It is now clear from (fTUp. ([Tlj) and that 

-i| < Ce~ Hu , (70) 
|. 9 -1| <Ce~ Hu , (71) 
|.g - 1 + Ar 2 /3| <Ce- Hu . (72) 

In particular, (|2"2"|) implies that 

m(u) < Ce~ Hu , 

and so the final Bondi mass M\ vanishes. Finally, geodesic completeness is easily obtained from (fT0")) - (fT2")) . 
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Appendix: the Einstein equations 

For the metric ([T]) we have the nonvanishing 
• (inverse) metric components: 



9uu = -gg , g ur = g ru = -g 

rr 9_ ru ur 

y •> y y •> 
g g 



= r ! 9<p<p = r sin 6 — sin Ogee 
1 1 Q 6B 



r 2 sin 2 ( 



sin 



Christoffel symbols = \g» v {d a g fiv + dpg va - d v g a0 ): 



"pU 

uu 

g 



d 19 



la/) — — , 1 ,„,„ — — Sin tfl (,(, 



9 



r r 1 9 i ~\ 9 



d Id. 
d^ 9 ~2d-r i99) 



T - = Tgd-r {99h T - = -gd-r 9 > 



T r 6e = - 9 -r, VI = -^sin 2 fl = sin 2 0rL 



sin o cos ( 



g 

r e - I - - — 



sin 



1 a <P<? ' 



Ricci tensor components R a p = d^T^ - d a T^ + T^F^ - T^T^ a : 



pr /pr 

uu \ rr 



Br 



r v ) 

(pr ) 



pr /p« 

ru \ uu 



- ru/ 



1?„ 



I?, 



1 / 9 2 



dud'. 



-(gg) 



2g \drdu 
1 / d . _ d , _ N n ~9g\ 

— -^-(gg) + g^-(gg) ~ 2 5 — 

r<? \ cm or aw J 



9d^(9~g)-2g- rdu: 



g + 



2r 



or du 



~dg d 
9 d-rd-r {99) 



d V r 

u r L ur 



8 r r 



Kr (n 



1 



o 2 



(99) 



o 2 



2g \dr 2 

~d r (rg r 

21dg 
r g dr ' 

+ 8 r T 



-^-N + i 



a P« 



-5 



dudr 
Y"P ) + T r (V 

(pr } 1 rr \ 



(pr) 
1 



<9r dr 



(gg) 



ipr; 



1 Sr 1 9r 



2 ^5% 
8u dr 

(pr (pr 



i_d_ 

rg dr 



(gg) 



d r r 

u r L a 



d e n 



deTl 



1 96\ l ru ' 1 uu) + 1 60 1 rr ~ 1 v 9 l V 8 



pu /pit I pr ) i pr pr _ p# pV 
(^(^p v- 1 - uu ' ur) ' (p(p rr 



sin 0i?ee 
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The Einstein field equations Q then have the following nontrivial components: 

1 f~d 2 , „_ d 2 \ 1 / <9<?<9<? „dg d , 
Y g \9^a{99) ~ 29g^9) + W , {*9frfo ~ 9^(99) 

r ou \g J rg or 



2g\dr i{99) 2 dudr g ) 2g 2 [dr dr^ 2 ' du dr 

1 d 

H — = k (d u <t>) (d r <j>) - Kg 

rg or 



(73) 



(74) 



2 -^ = K {d r(t >)\ (75) 
r g or 

^-(rg)=g(l-Ar 2 ) . (76) 



The wave equation reads 



- (a„ - r;j {d r <j>) + 9 -(d r - r; r ) (d r </>) - 4 Wr<« - 4 = o 

g g r z r z 



1 

r 



9 g 9 
du 2 <9r 



dr 2 \ dr J \ dr 



(77) 



References 



[And05] Michael T. Anderson, Existence and stability of even dimensional asymptotically de Sitter 
spaces, Annales Henri Poincare 6 (2005), 801-820, |a"rXiv:gr-qc/0408072| 

[Bey09a] Florian Beyer, Non-genericity of the Nariai solutions: I. Asymptotics and spatially homogeneous 
perturbations, Class. Quantum Grav. 26 (2009), 235015, arXiv:0902.253U 

[Bey09b] Florian Beyer, Non-genericity of the Nariai solutions: II. Investigations within the Gowdy class, 
Class. Quantum Grav. 26 (2009), 235016. larXiv:0902. 25321 

[Bey09c] Florian Beyer, The cosmic no-hair conjecture: a study of the Nariai solutions, Proceedings of the 
Twelfth Marcel Grossmann Meeting on General Relativity (R. Ruffini T. Damour, R.T. Jantzen, 
ed.), 2009. larXiv:1012.0056l 

[Bra97] Patrick Brady, Chris Chambers, William Krivan, and Pablo Laguna, Telling Tails in the Pres- 
ence of a Cosmological Constant, Phys. Rev. D 55 (1997), 7538-7545, arXiv:gr-qc/9611056 

[CAN12] Joao L. Costa, Artur Alho, and Jose Natario, Spherical linear waves in de Sitter spacetime, J. 
Math. Phys. 53 (2012), 024001. [arXiv:1107.0802l 

[Cho93] Matthew Choptuik, Universality and scaling in gravitational collapse of a massless scalar field, 
Phys. Rev. Lett. 70 (1993), 9-12. 

[Chr86] Demetrios Christodoulou, The problem of a self- gravitating scalar field, Commun. Math. Phys. 
105 (1986), 337-361. 

[Chr91] , The formation of black holes and singularities in spherically symmetric gravitational 

collapse, Commun. Pure Appl. Math. 44 (1991), 339-373. 

[Chr09] , The formation of black holes in general relativity, EMS Monographs in Mathematics, 

European Mathematical Society, 2009, larXiv:0805.3880l 

[Fri86] Helmut Friedrich, On the existence of n-geodesically complete or future complete solutions of 
Einstein's field equations with smooth asymptotic structure, Commun. Math. Phys. 107 (1986), 
587-609. 



22 



[GM07] Carsten Gundlach and Jose Martm-Garcfa, Critical Phenomena in Gravitational Collapse, Liv- 
ing Rev. Relativity 10 (2007), http://www.livingreviews.org/lrr-2007-5 

[LK11] Christian Lubbe and Juan Valiente Kroon, A conformal approach for the analysis of the non- 
linear stability of pure radiation cosmologies, arXiv:1111.4691 

[MN08] Hideki Maeda and Masato Nozawa, Generalized Misner- Sharp quasi-local mass in Einstein- 
Gauss-Bonnet gravity, Phys. Rev. D77 (2008), 064031. larXiv:0709.1199l 



[Nak95] Ken-ichi Nakao, On a quasilocal energy outside the cosmological horizon, arXiv:gr-qc/9507022 



[Ren04] Alan D. Rendall, Asymptotics of solutions of the Einstein equations with positive cosmological 
constant, Annales Henri Poincare 5 (2004), 1041-1064, |arXiv:gr-qc/0312020 

[Rin08] Hans Ringstrom, Future stability of the Einstein-non-linear scalar field system, Invent. Math. 
173 (2008), 123-208. 

[RS09] Igor Rodnianski and Jared Speck, The stability of the irrotational Euler- Einstein system with 
a positive cosmological constant, arXiv:0911.5501 

[Spell] Jared Speck, The nonlinear future- stability of the FLRW family of solutions to the Euler- 
Einstein system with a positive cosmological constant, accepted for publication in Selecta Math- 
ematica, 78 pages, arXiv:1102.1501. 

[TNN05] Sophonie Blaise Tchapnda N. and Norbert Noutchegueme, The surface- symmetric Einstein- 
Vlasov system with cosmological constant, Math. Proc. Cambridge Phil. Soc. 138 (2005), 541- 
553, |arXiv:gr-qc/0304098| 

[TNR03] Sophonie Blaise Tchapnda N. and Alan D. Rendall, Global existence and asymptotic behavior 
in the future for the Einstein- Vlasov system with positive cosmological constant, Class. Quant. 
Grav. 20 (2003), 3037-3049, |arXiv:gr-qc/0305059| 

[Wal83] Robert M. Wald, Asymptotic behavior of homogeneous cosmological models in the presence of 
a positive cosmological constant, Phys. Rev. D28 (1983), 2118-2120. 



23 



